Measurement of the Spectroscopy of Orbitally Excited B Mesons 

with the L3 detector 
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We measure the mass, decay width and production rate of orbitally excited B mesons in 1.25 
million hadronic Z decays registered by the L3 detector in 1994 and 1995. B meson candidates are 
inclusively reconstructed and combined with charged pions produced at the event primary vertex. 
An excess of events above the expected background is observed in the B7r mass spectrum near 
5.7 GeV. These events are interpreted as resulting from the decay B** — > B'*'-7r, where B** denotes 
a mixture of L = 1 B meson spin states. The masses and decay widths of the B2 (j q = 3/2) and 
BJ (j q = 1/2) resonances and the relative production rate for the combination of all spin states are 
extracted from a fit to the mass spectrum. 



I. INTRODUCTION 

Detailed understanding of the resonant structure of orbitally excited B mesons provides important information 
regarding the underlying theory. A symmetry (Heavy Quark Symmetry) arises from the fact that the mass of the b 
quark is large relative to Aqcd- In this approximation, the spin of the heavy quark (sq) is conserved independently 
of the total angular momentum (J q = s q + 1) of the light quark. Excitation energy levels are thus degenerate doublets 
in total spin and can be expressed in terms of the spin-parity of the meson J p and the total spin of the light quark 
j q . Corrections to this symmetry are a series expansion in powers of 1/toq [pL calculable in Heavy Quark Effective 
Theory (HQET). 

The L = mesons, for which j q — 1/2, have two possible spin states: a pseudo-scalar P (J p — 0~) and a vector V 
(J p = If the spin of the heavy quark is conserved independently, the relative production rate of these states is 
V/ (V + P) = 0.75.Q Recent measurements of this rate for the B system fls|-|ll"| agree well with this ratio. 

In the case of orbitally excited L = 1 mesons, two sets of degenerate doublets are expected: one corresponding to 
j q = 1/2 and the other to j q = 3/2. Their relative production rates follow from spin state counting (2J + 1 states). 
Rules for the decay of these states to the IS states are determined by spin-parity conservation For the dominant 
two-body decays, the j q = 1/2 states can decay via an L = transition (S-wave) and their decay widths are expected 
to be broad in comparison to those of the j q = 3/2 states which must decay via an L = 2 transition (D-wave). Table | 
presents the nomenclature of the various spin states for L = 1 B mesons containing either a u or d quark, with the 
predicted production rates and two-body decay modes. 

Several models, based on HQET and on the charmed L = 1 meson data, have made predictions for the masses and 
widths of orbitally excited B mesons. Some of these models |2|^,|^] place the average mass of the j q = 3/2 states above 
that of the j q = 1/2 states, while others [^,f7j predict the opposite ("spin-orbit inversion"). 



TABLE I. Spin states of the L = 1 mesons with their predicted production rates and decay modes. 



Name 


k 


J p 


Production 


Decay mode 


Transition 


BS 


1/2 


0+ 


1/12 


BS -> Btt 


S-wave 


BI 


1/2 


1+ 


3/12 


B* -» B*7T 


S-wave 


Bi 


3/2 


1+ 


3/12 


Bl -> B*7T 


D-wave 


B^ 


3/2 


2+ 


5/12 


B^ -» B*7T,B7T 


D-wave 



* Representing the L3 collaboration. 

1 Corrections due to the decay of higher excited states are predicted to be small. 
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Recent analyses at LEP combining a charged pion produced at the primary event vertex with an inclusively 
reconstructed B meson jfl]--[lj| have measured an average mass of Mb*» = 5700 — 5730 MeV, where B** indicates a 
mixture of all L = 1 spin states. An analysis [ fl4| combining a primary charged pion with a fully reconstructed B 
meson, measures Mb* = (5739 l^stat) ^(syst)) MeV by performing a fit to the mass spectrum which fixes the 
mass differences, widths and relative rates of all spin states according to the predictions of Eichten, et ai. 

The analysis presented here |jj is based on the combination of primary charged pions with inclusively reconstructed 
B mesons. Several new analysis techniques make it possible to improve on the resolution of the B-7T mass spectrum 
and to unfold this resolution from the signal components. As a result, measurements are obtained for masses and 
widths of D-wave Bj decays and of S-wave BJ; decays. 



II. EVENT SELECTION 
A. Selection of Z — > bb decays 

The analysis is performed on data collected by the L3 detector Jl5| in 1994 and 1995, corresponding to an integrated 
luminosity of 90 pb -1 with LEP operating at the Z mass. Hadronic Z decays are selected Jll| which have an event 
thrust direction satisfying | cos6*| < 0.74, where 8 is the polar angle. The events are also required to contain an event 
primary vertex reconstructed in three dimensions, at least two calorimetric jets, each with energy greater than 10 GeV, 
and to pass stringent detector quality criteria for the vertexing, tracking and calorimetry. A total of 1, 248, 350 events 
are selected. A cut on a Z — > bb event discriminant based on track DCA significances Jl7) yields a 6-enriched sample 
of 176, 980 events. 

To study the content of the selected data, a sample of 6 million hadronic Z decays have been generated with 
JETSET 7.4 (l]|, and passed through a GEANT based Jl9| simulation of the L3 detector. From this sample, the 
Z — > bb event purity is determined to be ir^ — 0.828. 



B. Selection of B" — > B^7r decays 

Secondary decay vertices and primary event vertices are reconstructed in three dimensions by an iterative procedure 
such that a track can be a constituent of no more than one of the vertices. A calorimetric jet is selected as a B candidate 
if it is one of the two most energetic jets in the event, if a secondary decay vertex has been reconstructed from tracks 
associated with that jet, and if the decay length of that vertex with respect to the event primary vertex is greater 
than 3cr, where a is the estimated error of the measurement. 

The decay of a B** to a B^*' meson and a pion is carried out via a strong interaction and thus occurs at the primary 
event vertex. In addition, the predicted masses for the L = 1 states correspond to relatively small Q values, so that 
the decay pion (tt**) direction is forward with respect to the B meson direction. We take advantage of these decay 
kinematics by requiring that, for each B meson candidate, there is at least one track which is a constituent of the 
event primary vertex and which is located within 90 degrees of the jet axis. A total of 60, 205 track-jet pairs satisfy 
these criteria. 

To decrease background, typically due to charged fragmentation particles, only the track with the largest component 
of momentum in the direction of the jet is selected. This choice has been found |^0|,[l4| to improve the purity of the 
signal. The track is further required to have a transverse momentum with respect to the jet axis larger than 100 MeV, 
to reduce background due to charged pions from D* — > Dir decays. These selection criteria are satisfied by 48, 022 Bit 
pairs with a b hadron purity of ttq = 0.942. 
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1. B meson direction reconstruction 

The direction of the B candidate is estimated by taking a weighted average in the 9 (polar) and <f> (azimuthal) 
coordinates of directions defined by the vertices and by particles with a high rapidity relative to the jet axis. A 
numerical error-propagation method pi| makes it possible to obtain accurate estimates for the uncertainty of the 
angular coordinates measured from vertex pairs. These errors, as well as the error for the decay length measurement 
used in the secondary vertex selection, are calculated for each pair of vertices from the associated error matrices. 

Particles coming from the decay of b hadrons produced in Z decays have a characteristically high rapidity relative 
to the original direction of the hadron when compared to that of particles coming from fragmentation. A cut on the 
particle rapidity distribution is thus a powerful tool for selecting the B meson decay constituents || . A second estimate 
for the direction of the B is obtained by summing the momenta of all charged and neutral particles (excluding the 
7r** candidate) with rapidity y > 1.6 relative to the original jet axis. Estimates for the uncertainty of the coordinates 
obtained by this method are determined from simulated B meson decays as an average value for all events. The final 
B direction coordinates are taken as the error-weighted averages of these two sets of coordinates. 

The resolution for each coordinate is parametrized by a two-Gaussian fit to the difference between the reconstructed 
and generated values. For 8 7 the two widths are <j\ = 18 mrad and 02 — 34 mrad with 68% of the B mesons in the first 
Gaussian. For (/>, the two widths are a\ = 12 mrad and cr 2 = 34 mrad with 62% of the B mesons in the first Gaussian. 



2. B meson energy reconstruction 

The energy of the B meson candidate is estimated by taking advantage of the known center of mass energy at LEP 
to constrain the measured value. The energy of the B meson from this method [l^| can be expressed as 

£ B =^ 2 + M B-^ coil , (1) 
2M Z V ' 

where Mz is the mass of the Z boson and M reC oii is the mass of all particles in the event other than the B. To 
determine M rGCO ii, the energy and momenta of all particles in the event with rapidity y < 1.6, including the it** 
candidate (regardless of its rapidity), are summed and M r 2 ecoil = Ey <16 — p y<16 . Fitting the difference between 
reconstructed and generated values for the B meson energy with an asymmetric Gaussian yields a maximum width 
of 2.8 GeV. 

III. ANALYSIS OF THE Bvr MASS SPECTRUM 

The combined B7r mass is defined as 



Mbtt = \JMl + ml + 2E B E„ - 2p B p^cosa , (2) 

where M B and m T are set to 5279 MeV and 139.6 MeV, respectively, and a is the measured angle between the B 
meson and the ir** candidate. The data mass spectrum is shown in Figure |l|.a together with the expected Monte 
Carlo background. 

A. Background function 

The background distribution is estimated from the Monte Carlo data sample, excluding B** — > B^tt decays, and 
fitted with a six-parameter threshold function given by 

p x x (x — p 2 ) P3 x e ( - PiX( - x ~ P2 ^ +P5X( - x ~ P2 ^ 2+P6><( - x ~ P2 ^^ . (3) 
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Parameters P2 through pe are fixed to the shape of the simulated background, while the overall normalization factor 
Pi is allowed to float freely in order to obtain a correct estimate of the contribution of the background to the statistical 
error of the signal. 

B. Signal function 

To examine the underlying structure of the signal, it is necessary to unfold effects due to detector resolution. The 
7r** candidates are expected to have typical momenta of a few GeV. In this range, the single track momentum 
resolution is no more than a few percent with an angular resolution better than 2 mrad. The dominant sources of 
uncertainty for the mass measurement are thus the B meson angular and energy resolutions. Monte Carlo studies 
confirm that these two components are dominant and roughly equal in magnitude. This analysis thus concentrates on 
unfolding the effects of these components by parametrizing and removing their contribution to the mass resolution. 

1. Signal resolution and efficiency 

The dependence of the B7r mass resolution and selection efficiency on Q value is studied by generating signal events 
at several different values of B** mass and Breit-Wigner width. The simulated events are passed through the same 
event reconstruction and selection as the data. The resulting B7r mass distributions are each fitted with a Breit-Wigner 
function convoluted with a Gaussian resolution (Voigt function) and the detector resolution is extracted by fixing the 
Breit-Wigner width to the generated value. 

The Gaussian width is found to increase linearly from 20 MeV to 60 MeV in the B** mass range 5.6 — 5.8 GeV. This 
increase with Q value is mainly due to the angular component of the uncertainty, which increases as a function of the 
opening angle a. The resolution is parametrized as a linear function of the B** mass from a fit to the extracted widths. 
Similarly, the selection efficiency is found to increase slightly with Q value and the dependence is parametrized with 
a linear function. 

Agreement between data and Monte Carlo for the B meson energy and angular resolution is confirmed by analyzing 
B* — > B7 decays selected from the same sample of B mesons. The photon selection for this test is the same as that 
described in reference ||. A B* meson decays electromagnetically and hence has a negligible decay width compared 
to the detector resolution. As in the case of the B-7T mass resolution, the B meson energy and angular resolution are 
the dominant components of the reconstructed B7 mass resolution. Fits to the Mb 7 — Mb spectra are performed 
with the combination of a Gaussian signal and the background function described above. For the Monte Carlo, the 
Gaussian mean value is found to be M Bj - Mb = (46.5 ± 0.6(stat)) MeV with a width of a = (11.1 ± 0.7(stat)) MeV. 
The input generator mass difference is 46.0 GeV. For the data, the Gaussian mean value is found to be Mb 7 — Mb = 
(45.1±0.6(stat)) MeV with a width of a = (10.7±0.6(stat)) MeV. Good agreement between the widths of the data and 
Monte Carlo signals provides confidence that the B energy and angular resolution are well understood and simulated. 

2. Combined signal 

According to spin-parity rules, five mass resonances are expected, corresponding to five possible B** decay modes: 
Bj — > B7r, Bj — > B*7r, Bi — » B*7r, B{ — » B*7r and Bq — > Bit. No attempt is made to tag subsequent B* — > B7 decays, 
as the efficiency for selecting the soft photon is relatively low. As a result, the effective B7r mass for a decay to a B* 
meson is shifted down by the 46 MeV B* — B mass difference. 
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FIG. 1. a) Mass spectrum for selected B7r pairs. The dots are data and the shaded histogram represents the expected 
background from Monte Carlo normalized to the sideband region 6.0 — 6.6 GeV; b) fit to the data Bit mass distribution with 
the five-peak signal function and the background function described in the text. 



The five resonances are fitted with five Voigt functions, with the relative production fractions determined by spin 
counting rules. The Gaussian convolutions to the widths are determined by the resolution function. Additional 
physical constraints are applied to the mass differences and relative widths in order to obtain the most information 
possible from the data sample. 

Predictions for the mass differences Mb* — M-q 1 and Mb* — Mb* depend on several factors, including the b and c 
quark masses and, in some cases, input from experimental data of the D meson system. The values are predicted to 
be roughly equal and in the range 5 — 20 MeV We constrain both of the mass differences to 12 MeV. 

Predictions for the Breit-Wigner widths of the j q = 3/2 are extrapolated from measurements in the D meson 
system Q and are expected to be roughly equal and about 20 — 25 MeV. No precise predictions exist for the 
j q = 1/2 states as there are no corresponding measurements in the D system. In general, however, they are also 
expected to be roughly equal, although broader than those of the j q — 3/2 states. We constrain T-q 1 = Tb* and 
Tb* = Tb* , but allow the widths of the and B^ to float freely in the fit. 



C. Fit results 



Monte Carlo events for each of the expected B** decays are generated and passed through the simulation and 
reconstruction programs and the B7r event selection. The resulting mass spectra are combined with background and 
fitted with the signal and background functions under the constraints described above. Mass values and decay widths 
for the BJ and B| resonances and the overall normalization are extracted from the fit and found to agree well with 
the generated values. All differences lie within the statistical error and have no systematic trend. 

The data Hit mass spectrum is fitted with the combined signal and background functions, allowing the normalization 
parameters to float freely. The resulting fit, shown in Figure [j], has a x 2 of 39 for 74 degrees of freedom. A total of 2652 
events occupy the signal region corresponding to a relative B** d production rate of cr(B** d )/cr(B u ,d) = 0.39±0.05(stat). 
The mass and width of the B^ are found to be M B * = (5770 ± 6(stat)) MeV and r B * = (21 ± 24(stat)) MeV and the 
mass and width of the B{ are found to be M B * = (5675 ± 12(stat)) MeV and r B * = (75 ± 28(stat)) MeV. 
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D. Systematic uncertainty 



Sources of systematic uncertainty and their estimated contributions to the errors of the measured values are summa- 
rized in Table ||. The b hadron purity of the sample is varied from 91% to 96%. The fraction of b quarks hadronizing 
to B Ui d mesons is taken to be 79% and is varied between 74% and 83% in accordance with the recommendations of 
the LEP B Oscillation Working Group |22[ , These variations effect only the overall B** production fraction. 

Systematic effects due to background modelling are studied by varying the shape parameters of the background 
function and by performing the fit with other background functions to study the effect on the measured values. 
Contributions to the error due to modelling of the signal are estimated for the mass and width constraints: the 
Mb« — Mb x and Mb* — Mb* mass differences are varied in the range 6 — 18 MeV and the Tb ± /Tb*. and Tb* /Lb* ratios 
are varied between 0.8 and 1. Effects due to uncertainty in the resolution and efficiency functions are estimated by 
varying the slopes and offsets of the linear parametrizations. 

Three-body decays of the type B2 — > Jiirir have been generated and passed through the simulation and reconstruction 
programs and the B7r event selection. B7r pairs, for which only one of the pions is tagged, are studied as a possible 
source of resonant background. The resulting reflection is found to contribute insignificantly to the background 
in regions of small Q value. Similarly generated B** — > BK decays, where the K is mistaken for a ir are found 
to contribute only slightly to the low Q value region and their effects are included in the background modelling 
uncertainty contribution. 



IV. CONCLUSION 

We measure for the first time the masses and decay widths of the B2 {j q = 3/2) and B| (j q = 1/2) mesons. From 
a constrained fit to the B7r mass spectrum, we find 

Mb* = (5770 ± 6(stat) ± 4(syst)) MeV 

T B ; = (23 ± 26(stat) ± 15(syst)) MeV 

AIb* = (5675 ± 12(stat) ± 4(syst)) MeV 

F B * = (76 ± 28(stat) ± 15(syst)) MeV . 

The relative B** d production rate, including all L = 1 spin states, is measured to be 

Br(b -> B** -> bMtt) 

— „ „ ' r = 0.39 ± 0.05(stat) ± 0.06(syst) 

Br(b -> B u , d ) 

where isospin symmetry is employed to account for decays to neutral pions. 



TABLE II. Sources of systematic uncertainty and their estimated contributions to the errors of the measured values. 



Sources 






Mb- 


r B * 


r 


b purity 










±0.02 


B Uj d fraction 










±0.03 


background 


±2 


±9 


±3 


±9 


±0.05 


M constraints 


±3 


±7 


±3 


±7 


< 0.01 


r constraints 


< 1 


±2 


< 1 


±2 


< 0.01 


resolution 


±2 


±9 


±1 


±9 


±0.01 


efficiency 


< 1 


±1 


< 1 


±1 


±0.01 


Total 


±4 


±15 


±4 


±15 


±0.06 
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